Domain Formation and Orbital Ordering Transition in a Doped Jahn- Teller Insulator 
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The ground state of a double-exchange model for orbitally degenerate e g electrons with Jahn- Teller 
lattice coupling and weak disorder is found to be spatially inhomogeneous near half filling. Using a 
real space Monte-Carlo method we show that doping the half-filled orbitally ordered insulator leads 
to the appearance of hole-rich disordered regions in an orbitally ordered environment. The doping 
driven orbital order to disorder transition is accompanied by the emergence of metallic behavior. We 
present results on transport and optical properties along with spatial patterns for lattice distortions 
and charge densities, providing a basis for an overall understanding of the low doping phase diagram 
of Lai-zCa^MnOs. 
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Hole-doped perovskite manganitcs, for example 
Lai-^Ca^MnOs (LCMO), are well known for their colos- 
sal magnetoresistance (CMR) effect 0, The 'opti- 
mally doped' CMR compounds with x ~ 0.3 have been 
the focus of numerous theoretical studies and are qualita- 
tively understood in terms of the interplay of the double 
exchange mechanism, electron-phonon interactions, and 
disorder in a single electronic band Q. 0- Less at- 
tention has been given to the low-doping regime, where 
the two-band character of manganites is crucial and in 
addition to charge, spin, and lattice variables, the or- 
bital degrees of freedom and their ordering become im- 
portant. The undoped compounds are orbitally ordered 
(00), A-type antiferromagnetic insulators with large 
Jahn- Teller (JT) distortions of the Mn06 octahedra p|, 
which lift the degeneracy of the two Mn-e 9 levels. Elec- 
tronic and cooperative lattice effects lead to a staggered 
ordering of the lattice distortions and the orbital occu- 
pancies. Upon doping, the antiferromagnetic insulator 
evolves into a ferromagnetic insulator, with weakened 
orbital order, and eventually undergoes a transition to 
an orbitally disordered ferromagnetic metal (OD-FM-M) 
0- Despite the achieved progress towards an under- 
standing of mag netic and orbital ordering in the undoped 
compounds [T(j, efforts for analyzing the doping driven 
transition from the orbitally ordered insulating to the or- 
bitally disordered metallic phase have remained limited. 
A simple view of this transition rests on an entirely clas- 
sical picture in terms of random fields introduced by the 
doped holes [Tl| . 

The critical hole doping xo d for the loss of orbital or- 
der is close to the doping ximt for the insulator- metal 
transition (IMT) in LCMO, where x IM t ~ 0.22 In 
lower bandwidth materials like Pri-a-Ca^MnOa (PCMO) 
the insulating phase persists to even larger hole concen- 
trations 01 ■ NMR and neutron scattering experiments 
suggest, that the doping regime x < ximt is spatially 
inhomogeneous in LCMO with coexisting 'hole poor' or- 
bitally ordered and 'hole rich' orbitally disordered regions 
|l4j . and the observation of confined 'spin waves' confirms 



the existence of magnetic clusters on the nanoscale |15( . 
It has remained unclear how the JT insulator evolves 
from the homogeneous OO state at x = to the homo- 
geneous OD metal at optimal doping through an inter- 
mediate inhomogeneous state. 

In this Letter we present results on a two-band double- 
exchange model with electron-lattice coupling and disor- 
der in two dimensions (2D) using a real-space technique. 
We provide a description for the doping driven loss of 
orbital order and the detailed doping vs. temperature 
phase diagram at intermediate electron-lattice coupling, 
appropriate to LCMO. Results for charge transport and 
spectral properties are presented, which characterize the 
metal-insulator transitions. Real-space structures for the 
inhomogeneous state of the hole doped insulator allow us 
to follow the emergence of orbitally disordered domain 
walls at low doping and their eventual percolation. 

Specifically, we consider a two-band model for itiner- 
ant e g electrons coupled to JT lattice distortions and to 
localized 5 = 3/2 ti g spins in the presence of substitu- 
tional disorder, described by the Hamiltonian: 

a/3 

{ij)a i i 

+J S s * • S J + di n + f E IQ«I 2 - W 

(ij) » » 

The magnetic properties arise from the competition be- 
tween the Hund's rule coupling Jy driven double ex- 
change and the antiferromagnetic superexchange Js be- 
tween the tig core spins S^. 

In Eq. (1), c and arc annihilation and creation oper- 
ators for e g electrons and a, (3 are summed over the two 
Mn-e 9 orbitals d x 2_ y 2 and d 3z 2_,.2, which arc labelled (a) 
and (b) in what follows. t"P are the hopping matrix ele- 
ments between e g orbitals on nearest-neighbor sites and 
have the cubic perovskite specific form: t% a = ty a = t, 

= th y = */3, tf = t b x a = -t/V3, tf = t b y a = t/y/3 
[l|, where x and y denote the spatial directions on a 
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square lattice. The disorder is modelled by random on- 
site potentials e,-, with equally probable values ±A. The 



e„-electron spin is cr£ 



"iaaKa'dcto', where F M 
are the Pauli matrices. A denotes the strength of the JT 
coupling between the distortion Qi = (Qi X , Qi Z ) and the 
orbital pseudospin if = YTf Aa^ap ^ Q- A' controls 
the lattice stiffness, and [i is the chemical potential. 

We set t = 1 as the reference energy scale. In the man- 
ganites Jh 3> 1, and we adopt the frequently used limit 
Jh — * oo, which retains the essential physics pj. We use 
J s = 0.05 throughout, which is estimated from the Neel 
temperature for CaMnC>3, where antifcrromagnctism is 
purely superexchange driven. The parameters A and A 
will be selectively explored. The spins are assumed to 
be classical unit vectors, |S, | = 1; quantum effects in the 
lattice variables are not considered, and the stiffness is 
set to K = 1 . In the limit Jh — > oo the spin of the e g 
electrons is tied to the orientation of the local core spin 
leading to a two-orbital 'spinless' fermion model with core 
spin configuration dependent hopping amplitudes 0. 

Replacing a fraction x of rare earth ions with 2 + 
cations in the parent manganites affects the mean A- 
site ionic radius Ta(x) as well as its variance <ja{x)- The 
varying ta modifies the electronic hopping amplitude, 
and hence the X/t ratio, while a a controls the disorder 
strength A. In most of what follows we set A = 1.6, 
which reproduces the transport gap ~ 0.4cV (if we as- 
sume t ~ 0.2eV) in LaMnOj, estimated from the acti- 
vated resistivity behavior [l(j. We set A = 0.4 as a 
typical value for weak disorder, and explore the doping 
and temperature dependence. Naturally the amount of 
disorder depends on the doping level, but this variation 
is not addressed here. 

The model defined in Eq. (1) has been studied ear- 
lier using mean-field methods, as well as exact diagonal- 
ization (ED) based Monte-Carlo (MC) simulations [l7j . 
While the mean-field approximation excludes by con- 
struction the possible existence of inhomogeneous phases, 
the accessible system sizes within ED-MC are too small 
(~ 100 sites) to explore spatial clustering effects, the or- 
bital order to disorder transition, or the IMT itself. Here 
we use the travelling cluster approximation (TCA), which 
readily allows access to systems of ~ 1000 sites to anneal 
the classical spin and lattice variables. ED of the full 
fermionic Hamiltonian is used only for computing the 
electronic quantities in the TCA-generated classical con- 
figurations. This method has been benchmarked before 
[l8j and applied to a one-band version of Eq. (1) ||. 

Fig. 1 summarizes our results for the OO-OD transi- 
tion, the core-spin magnetism, and the doping and tem- 
perature driven IMT. Panel (a) shows the x — T phase 
diagram for A = 1.6 and A = 0.4. Squares mark the fer- 
romagnetic (FM) to paramagnetic (PM) crossover and 
circles denote the orbital ordering transition tempera- 
tures Too as inferred from the temperature dependence 
of the q = (it, it) = qo component of the lattice structure 




factor, D(, (q) = ^ Eij(QrQ 
in panel (c) . Here and below 



j/av c !q '( ri r ^ shown 
denotes the combined 



FIG. 1: (Color online) (a) x — T phase diagram at A = 1.6. 
(b) x — X phase diagram at T = 0.02. The disorder strength is 
A = 0.4. The shaded region in (b) indicates phase separation. 
OO (OD) refers to an orbitally ordered (disordered) phase; 
FM (PM) denotes a ferromagnetic (paramagnetic) state and 
M (I) indicates metallic (insulating) character. T dependence 
of (c) the lattice structure factor Dq{(tt, n)) and (d) the den- 
sity of states at the chemical potential N(fj.) for different x. 



average over thermal equilibrium configurations and over 
the realizations of quenched disorder. The I-M boundary 
is obtained from the sign of the slope of the resistivity 
p(T) (sec Fig. 2(a)). At T = the system is ferro- 
magnetic at all doping levels, despite the presence of the 
antifcrromagnetic superexchange coupling Jg; the dop- 
ing driven OO-OD transition at xod ~ 0.22 occurs close 
below the IMT. For x < xod the system loses either 
ferromagnetic order (for x — > 0) or orbital order first 
(for x near xod) with increasing T, and for T > 0.1 all 
long range order is lost leading to an orbitally disordered 
paramagnetic insulator (OD-PM-I). 

A broader perspective for the OO-OD transition is ob- 
tained from the low temperature A — x phase diagram for 
A = 0.4, shown in Fig. 1(b). For our choice of Js the 
system is FM over the entire selected parameter range. In 
the clean limit, A = 0, the major feature is phase separa- 
tion (PS), as indicated by the grey shaded area, between 
the OO-FM-I at x = and the OD-FM-M for x > 0. PS 
is identified from the existence of a jump in the average 
hole density upon varying the chemical potential. In the 
presence of disorder the PS range is replaced by inho- 
mogeneous phases. We identify three distinctly different 
regimes with respect to A: (i) For A < 1 orbital ordering 
is absent and the system is metallic over the entire dop- 
ing range, (ii) When 1 < A < 1.6 the system sustains 
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orbital order even away from x — 0, and the orbital order 
and the insulating character is lost at a critical doping 
xod(X) ~ Ximt(X), which reflects that the insulating 
character is intimately related to the staggered ordering 
of the lattice distortions. (Hi) For A > 1.6 the system 
first loses orbital order and subsequently becomes metal- 
lic at higher doping, so ximt > xod- In fact, for A > 2 
the system does not become metallic even for x > 0.3. 
This is the strong coupling regime where single electrons 
can be 'self-trapped' and orbital ordering is not a pre- 
requisite for the insulating behavior. If we assume that 
X/t ~ 1.6 is appropriate to LCMO, the phase diagram 
in Fig. 1(a) agrees remarkably well with the experimen- 
tal results reported in Ref. 12]. The lower bandwidth 
materials would correspond to larger X/t, and indeed we 
find that the insulating character persists to much larger 
x for larger A, consistent with the experiments on PCMO 
[lflj . We remark that some aspects of the 3D manganite 
physics are not contained in our 2D calculation, e.g. we 
can not address the crossover from an A-type antiferro- 
magnet to a FM state; the present results should thus 
be compared to the physics of a single plane of the real 
materials only. 

Fig. 1(c) shows the temperature dependence of 
Dg(qo)j which is a measure of the staggered ordering 
tendency of the JT distortions, for different hole concen- 
trations. Since the lattice distortions are coupled to the 
orbital pseudospin, the same tendency is transferred to 
the analogously defined orbital structure factor D r (qo). 
Therefore Dq(c\q) serves as an indicator for the stag- 
gered ordering of both, the lattice distortions and the 
orbital pseudospin. Staggered order is found only in the 
Q x component of the distortions. 

In the undoped case the onset of orbital order is 
accompanied by the opening of a gap around the 
Fermi level in the density of states (DOS) N(u>) = 
N- 1 (En S & - e n{Q, S})) av , where e„{Q, S} denote the 
single particle eigenvalues for a configuration {Q,S} of 
the classical variables. Fig. 1(d) tracks the tempera- 
ture dependence of N(p) for the same choice of x values 
as in Fig. 1(c). N(p) vanishes at low temperatures for 
x < 0.10, indicative for insulating behavior but retains 
a small finite value for x > 0.10; the latter is in fact re- 
lated to a pseudogap structure in the DOS and suggests 
a possibly 'metallic' state (see Fig. 2). 

In order to discuss the implications for transport we 
compute the optical conductivity a(u) using the Kubo 
formula with the exact eigenstates [lj|- The resistiv- 
ity p is approximated by the inverse of cr(u) m i n ), where 
w m j n = 20t/N ~ 0.03t is the lowest reliable frequency 
scale for <j(uj) calculations on our N = 24 2 system. Fig. 
2(a) shows the temperature dependence of p for different 
x. For x < 0.2 there is a sharp rise in p(T) at the on- 
set temperature for orbital ordering Too- For x ~ 0.3 
and x ~ 0.4, however, there is a downturn in p(T) upon 
cooling. We correlate this behavior with the tempera- 
ture dependence of the magnetization m(T) defined via 
m 2 = ((N- 1 J2i S 2 ) atl > shown in Fi S- 2 ( b )- For a11 val " 
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FIG. 2: (Color online) (a) Resistivity p on a logarithmic scale 
in units of h/ne 2 and (b) magnetization m as a function of 
temperature for different x. The low temperature (c) optical 
conductivity cr(a>) and (d) density of states N(u>) for varying 
x. The results are for A = 1.6 and A = 0.4. A Lorentzian 
broadening of 0.02 is employed for the DOS calculations. 



ues of x a para- to ferromagnet transition is observed 
upon cooling, with the Curie temperature T c inferred 
from the inflection point in m(T). Note that T c does 
not increase with increasing x, contrary to the simple 
double exchange scenario |20|. Fig. 2(c) shows the opti- 
cal conductivity at low temperature for different x. For 
x < 0.2 the low- frequency spectral weight is strongly sup- 
pressed, which is evident in the effective carrier density 
n e ff(cu,x) = Jq a(uj')dio' at u = 1, shown in the inset. 
For x = 0.3, er(cj) remains finite at the lowest attainable 
u), but the response is non-Drude like, while at x — 0.4 
a more conventional optical response emerges. Neverthe- 
less, even at x = 0.4, n e // is strongly suppressed com- 
pared to the naive electron count. This is consistent with 
our observation in the disordered one-band model Q and 
appears as a generic feature of the interplay of disorder 
and electron-lattice coupling. Fig. 2(d) shows how the 
DOS evolves from the 'clean gap' at x = through the 
low x pseudogap to the high T c 'metallic' regime. The 
clean gap at x = 0, which or igin ates from the nesting in- 
stabilities at weak coupling |2l| and effective repulsions 
between self-trapped electrons at str ong coupling, is sta- 
ble in the presence of weak disorder |23 ■ From the com- 
bination of N(p), p(T), and <j(lo) we conclude that the 
electron system at T = has an TMT' near x = 0.25. 

For a more microscopic understanding of these re- 
sults, we have examined the doping evolution of the real 
space patterns for the local charge density rii and the 
local spatial correlations of the lattice variables Cq = 

I Ea Qi ' Qi+i5 f° r a representative disorder realization, 
where 8 is summed over the nearest neighbors of site i. 
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FIG. 3: Spatial patterns at x = 0.04, 0.11 and 0.19 (left, mid- 
dle and right). Top row: charge density m, grayscale covers 
the range from 0.2 (black) to 1 (white). Bottom row: local 
correlations of the lattice distortions Cq (see text). Grayscale 
from -1.5 (white) to (black). All results are for a specific 
disorder realization on a 40 x 40 lattice, T — 0.01, A = 1.6 
and A = 0.4 



The rii in the upper row in Fig. 3 show (with white 
marking the 'hole poor' and black the 'hole rich' regions) 
that the doped holes at low x lead to a strong density 
variation, although the holes are not 'site localized'. The 
hole positions are spatially correlated with a short range 
charge order pattern. The spatial pattern is filamentary, 
rather than 'puddle-like', with the linear structures con- 
necting up for x ~ xod- The Cq pattern is understood 
from the hole locations, and the above discussed supres- 
sion of .Dq (qo) arises from the loss of OO in the vicinity 
of the holes as well as from the presence of antiphase 
domains separated by the 'hole rich' domain walls. 

Although local Coulomb interactions were not explic- 
itly included in our model analysis their effects are never- 
theless partially captured. E.g. the large Jh avoids dou- 



ble occupancy of a single orbital and therefore acts like an 
intra-orbital Hubbard repulsion. The JT polaron binding 
energy, ~ X 2 /2K, has effects similar to an inter- orbital 
Hubbard repulsion since it prefers one e g electron per 
site |lOj . Therefore we do not expect qualitative changes 
in our phase diagrams if explicit electron-electron inter- 
actions were included. The Too scale however will be 
affected. The cooperative nature of the lattice distortion 
will also enhance Too &s we have checked. The criti- 
cal doping for the OO-OD transition, however, does not 
seem to be significantly affected by cooperative effects, 
and neither are the spatial patterns for x ~ xod H3 ■ 

Naturally the 'IMT' that we observe in our 2D model 
with disorder is to be understood as a crossover from an 
insulating phase to a weakly localized (WL) phase with a 
finite density of states at the Fermi level. The localization 
effects become apparent by studying progressively larger 
lattices, and we observe the expected slow growth of the 
computed resistivity In 3D the gapped insulator to WL 
crossover is expected to become a genuine IMT. 

In summary, our results on the 2D Jahn- Teller double 
exchange model reveal a doping driven transition from 
an orbitally ordered insulator to an orbitally disordered 
ferromagnetic metal in agreement with the experiments 
on Lai-aCazMnO,-?. In the orbitally disordered regime, 
the system undergoes a thermally-driven transition from 
a ferromagnetic metal to a paramagnetic insulator, char- 
acteristic of the CMR materials. The intermediate inho- 
mogeneous phase, with coexisting orbitally ordered and 
orbitally disordered regions, allows a natural interpreta- 
tion of the neutron scattering and NMR data in LCMO. 
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